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Abstract A key ingredient of quantum repeaters is en¬ 
tanglement distillation, i.e., the generation of high-fidelity 
entangled qubits from a larger set of pairs with lower 
fidelity. Here, we present entanglement distillation pro¬ 
tocols based on qubit couplings that originate from ex¬ 
change interaction. First, we make use of asymmetric bi¬ 
lateral two-qubit operations generated from anisotropic 
exchange interaction and show how to distill entangle¬ 
ment using two input pairs. We furthermore consider the 
case of three input pairs coupled through isotropic ex¬ 
change. Here, we characterize a set of protocols which 
are optimizing the tradeoff between the fidelity increase 
and the probability of a successful run. 


1 Introduction 

In a quantum communication (QC) network, the estab¬ 
lishment of long-distance entanglement is indispensable 
to fully harness the advantages offered by quantum infor¬ 
mation processing [Tj, e.g. perfectly secure long-distance 
quantum communication [2]. For the distribution of max¬ 
imally entangled states, one has to counteract decoher¬ 
ence processes due to the unavoidable interaction of en¬ 
tangled particles with their environment. A fundamen¬ 
tal component of a QC network are therefore quantum 
repeaters 13 i, which enable the successive creation 
of near-maximal entanglement between distant network 
nodes. Entanglement distillation, on the other hand, is a 
key part of quantum repeaters, and requires a function¬ 
ing quantum memory [3- 

Spins in solid-state environments, such as single elec¬ 
trons in semiconductor quantum dots (QDs) 0 or nitrogen- 
vacancy centers in diamond .fj, show long coherence 
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times (fj .s to ms) and offer flexible controlling mecha¬ 
nisms. However, the original protocols of entanglement 
distillation 0 1 are rather unpractical, e.g. for spin 
qubits in QDs mentioned above, since an efficient imple¬ 
mentation of the required controlled-NOT (cnot) gates 
is very demanding. In the case of Heisenberg exchange, 
it requires two two-qubit interaction pulses each leading 
to the yj SWAP gate, and additionally five single qubit ro¬ 
tations to construct a CNOT gate m- However, single¬ 
spin rotations take on the order of 100 ns d! and are 
therefore slow compared to exchange-based two qubit 
operations. The ^SWAP gate, e.g., has been succesfully 
implemented in less than 200 ps [T2] . 

This circumstance motivates the work presented in 
this contribution, namely a careful study of entangle¬ 
ment distillation protocols using only the typical inter¬ 
action between electrons in QDs, namely the exchange 
interaction mm- In the following, we first extend an 
earlier proposal based on isotropic Heisenberg exchange 
M , where the concept of asymmetric bilateral two-qubit 
operations for protocols using two input pairs was intro¬ 
duced, to the more general scenario of an anisotropic 
exchange interaction. 

Furthermore, we analyze entanglement distillation pro¬ 
tocols for exchange-coupled qubits that use three input 
pairs. Our method is based on an algebraic view of the 
occurring operations and we find protocols optimizing 
the tradeoff between the gain in fidelity and the proba¬ 
bility of a successful run. 

2 Preliminary remarks 

An orthonormal basis of the two-qubit Hilbert space is 
given by the maximally entangled Bell states, 

|^ ± ) = ^ (|tt) ± | U )) , (1) 

|^ ± ) = 2=du)±ut)), (2) 
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where the two spin eigenstates | f) = 10) (up) and 14-) = 

11) (down) define the computational basis. The overlap 
F of an arbitrary two-qubit quantum state p with the 
Bell state | ), i.e. 

F=(F-\p\F-), (3) 

is referred to as the fidelity of the state p in the following. 

Recurrence protocols work on two or more qubit pairs 
of low fidelity as input that are used to create a single 
qubit pair with higher fidelity as output. Thereby, only 
local unitary operations, measurements, and two-way 
communication of the measurement results via a clas¬ 
sical channel can be used. Having initially many copies 
of the low-fidelity pairs and running the distillation pro¬ 
tocol iteratively on the output pairs with higher fidelity, 
one can achieve fidelities arbitrarily close to F = 1 and 
thus, obtain a maximally entangled state. 

The original idea of entanglement distillation was in¬ 
troduced in Ref. [S) and will be referred to as the BBPSSW 
protocol in the following. Initially, the physical setup is 
such that the two communicating parties, commonly re¬ 
ferred to as Alice and Bob, have access to mixed two- 
qubit states pi of fidelity F = (\ p l \ ) < 1 that 

can originate from imperfect sources or noisy quantum 
channels. To apply the distillation protocol, the state pi 
first needs to be brought into the Bell-diagonal form 

Pf = F |<2>+)(<2> + | 

+ 1 -^ F ( | F + ) ( F + | + | | + | ) (<2>~ | ) . 

( 4 ) 

This can be achieved for an arbitrary two-qubit state by 
a so-called twirl operation Hus that retains the compo¬ 
nent of the rotationally invariant state | ), equalizes 

the components of the other three Bell states, and re¬ 
moves all off-diagonal elements. Thereby Alice and Bob 
have to implement a random bilateral rotation, i.e. they 
choose a random SU(2) rotation and apply it locally to 
each of the qubits, respectively. As an intermediate re¬ 
sult, a so-called Werner state W F [IS] is created, 

W f = F \&-)(&~\ 

+ (^+1 + |<z> + )( <? + I + |<2>“ ) (<z>~ I ), 

( 5 ) 

that can be brought into the form in Eq. Q by per¬ 
forming a unilateral rotation of tt about the y axis on 
the Bloch sphere of one of the two qubits, thereby inter¬ 
changing the | ) and | <P + ) components. 


3 Asymmetric entanglement distillation with 2 
pairs of spins 

3.1 General interactions between two spins 1/2 
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Figure 1 Entanglement distillation with the exchange in¬ 
teraction. (a) Two input pairs. The parties Alice and Bob 
share two imperfectly entangled input pairs p F , qubit labels 
are specified. Different exchange pulses generated by the in¬ 
teractions J A (t) (Alice) and J B (t) (Bob) implement an asym¬ 
metric bilateral two-qubit operation, (b) Three input pairs. 
Alice and Bob have control on interactions between next lo¬ 
cal neighbors. Protocols are implemented by iteration of local 
two-qubit operations. 


symmetric bilateral CNOT gate by an asymmetric bilat¬ 
eral operation, in which each local two-qubit operation 
between qubits i and j is generated from the (an)isotropic 
exchange interaction 

H a(t) = ( a x )a x ] + a y )(J y ) + > ( 6 ) 


where the crft (p = x,y,z) are the Pauli matrices de¬ 
scribing the ith qubit. The parameter £ quantifies the 
anisotropy of the interaction, e.g. for £ = 1 the Hamilto¬ 
nian in Eq. (6j) describes isotropic exchange interaction 
(see Sec. 
i£] 


3.21. The time evolution generated by Hij(t) 


Uij(a ) 


— F. 4 /o Hij (t ) 


fe~ iS f 0 0 0 \ 

0 e i5 rcos(f) sin (§) 0 

0 -ie^sind) e 4 ^ cos (f) 0 

0 0 0 e -4 T / 

( 7 ) 


Here, the matrix representation is in the product basis 
{|Tt) i Itl) > lit)! Ill)} and we set the initial time to 
zero. The time evolution is parametrized by the so-called 
pulse area a defined as 


a = 


t 


J 


d t 1 


( 8 ) 


For the distillation of partially entangled states that 
are produced by a source and subsequently transmitted 
to Alice and Bob, they can use two copies of the state 


We start our description with a distillation scheme sim- 1 Here, we set h = 1, and time-ordering in Eq. |t| is not 
ilar to the original BBPSSW protocol 0. We replace the necessary since [//,;,(t), = 0 for all t and t'. 
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PF [Eq. @L which can be produced from arbitrary two- 
qubit states with fidelity F as described in Sec. [2j The 
qubits at Alice’ site are labelled 1 and 3, and Bob pos¬ 
sesses qubits 2 and 4 (see Fig. [TJ. The four-qubit state 
p describing such a system is thus given by 


( 12 ) 

P = Pf 




(9) 


where p denotes the state of qubits i and j. After¬ 
wards, Alice and Bob each apply an exchange pulse be¬ 
tween their respective qubits, which is described by the 
unitary transformation 

U(a,/3) = U 13 (a) <g> U 2 i{0). (10) 


If we denote the exchange couplings in Alice’s and Bob’s 
spin register as J A (t) and J B (f), then the pulse areas a 
and /3 are given by 

dt'J A (t'), (11) 

P= r dt ( 12 ) 

Jo 

where t A / B are the respective pulse lengths. The cru¬ 
cial difference to other entanglement distillation proto¬ 
cols HM2 is that Alice and Bob are allowed to choose 
different pulse areas and thus, apply different bilateral 
two-qubit operations. It is exactly this asymmetric bilat¬ 
eral operation that makes entanglement distillation via 
exchange interaction of the form in Eq. § feasible at 
all if only two input pairs are used. The exchange pulses 
transform the four-qubit state p as 

p^U(a,p) P U(a,p)l (13) 



After this unitary transformation, the two parties con¬ 
tinue in the same way as in the original BBPSSW pro¬ 
tocol. Although we do not use any conditional quan¬ 
tum operations here, we still denote qubits 1 and 2 as 
control qubits, and qubits 3 and 4 as target qubits. Al¬ 
ice and Bob measure the target qubits in the compu¬ 
tational basis {|t) j |l)} and compare the measurement 
results afterwards using classical two-way communica¬ 
tion. If Alice and Bob obtain equal measurement results, 
i.e. either both spins are pointing up or both are point¬ 
ing down, they will keep the control qubits. Otherwise, 
the state is discarded. In case that the control qubits 
are kept, another unilateral rotation of n about the y 
axis on the Bloch sphere is applied to interchange again 
the | } and | <P + } components. As we derive below, in 

case of keeping the control pair, the fidelity of precisely 
this state can become larger than the initial fidelity F 


2 We can separate the time evolution of the four-particle 
system in Eq. (10 1 into the two-particle propagators [/ 13 (a) 
and U24(l3) because the Hamiltonians describing each ex¬ 
change interaction commute, i.e. [H\ B (t), H 24 (t')\ = 0 for all 
t and t', with Hij(t ) given in Eq. (| 6 |. 



Figure 2 Increased fidelity F ou t(F, a, —a) after a single dis¬ 
tillation step as a function of F, here shown for pulse areas 
a = 7t/4 and a = 7 r/ 2 . 


through the above transformation and measurement, de¬ 
pending on the applied exchange pulses a and (3. 

If we denote the postselected state of the control 
qubits by p', then the output fidelity F out (F, a, f3) = 
(| p' | ) is found to be 


F out (F,a,(3) 


v{F,ot,P) 

S(F,a,0y 


(14) 


with 

v(F, a , /?) = 3(4F — 1) cos(cr) cos(/3) 

+ 4(8F 2 + 2F-1) cos cos 

- (4 F - l) 2 sin(a) sin(/?) + 4F(4F + 1) + 7, (15) 


and 

S(F, a, (3) = 6(4 F — 1) cos(cr) cos(/3) 

— 2(4F — l) 2 sin(a) sin(/3) + 6(4F + 5) (16) 

A detailed analysis of the fidelity F out (F,a, (3) may be 
found in Ref. mi for the isotropic case £ = 1. An in¬ 
teresting property is found when Alice and Bob apply 
mutually inverse operations, i.e. a = —f3. In this case, 
the fidelity F out (F, a. —a) becomes independent of the 
anisotropy parameter £, 

/ i 1 

F ou t(F , a, —a) = - 

3 - 12E 2 

+ (F - 1)(4 F - 1) cos(2cr) - F(4F + 7) - 7 ^ 

Since in this case, a repulsive fixed point F m \ n = 1/2 
and an attractive fixed point F max = 1 of the map 
-Pout (F, a, —a) are found, it allows the distillation of max¬ 
imally entangled states. In the range 1/2 < F < 1, the 
maximum of F out (F , a,—a) is obtained for a = 7t/2. 
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3.2 Heisenberg exchange interaction 


Isotropic exchange interaction (£ = 1) is described by a 
Heisenberg Hamiltonian, i.e. 

' <*j, ( 18 ) 

where <x denotes the vector of Pauli matrices. Typically, 
electron spins in gate-defined quantum dots are coupled 
by such an interaction that can be used to implement 
universal quantum computation mmmmm- Tins 
case has been studied in Ref. HI and it was found that 
the highest gain in fidelity is given for pulse areas a = 
tt/2 = -/?, 


F out , opt (F) = -FLut (r, ^ + 27m, + 2mnJ 

_ 16F 2 + F + 1 
“ 8F 2 + 2F + 8 ’ 


(19) 


where the integers n and m can be chosen independently 
by Alice and Bob. In the optimal case, Alice applies the 
so-called V SWAP gate, 


U 


/1 0 0 0\ 

0 (1 — i)/2 (1 + i)/2 0 
0 (l+i)/2 (l-i)/2 0 ’ 

\0 0 0 1 / 


and Bob the inverse 


Vs wap gate, a/swap ^ 


( 20 ) 


3.3 XY interaction 


The Hamiltonian describing XY-type interaction is ob¬ 
tained for 5 = 0 and thus given by 

= \j{t) + <J y )(T y ) ) , ( 21 ) 


i.e. only the x and y components of the spins are cou¬ 
pled. This kind of interaction appears, e.g., in all-optical 
cavity-coupled QD electron spins .23] or superconduct¬ 
ing qubits |25j . For a pulse area a = — n, the Hamiltonian 
HfJ ( t ) generates e.g. the so-called iswAP gate, 


/I 0 0 0\ 
0 0*0 
0 * 0 0 
\0 0 0 1 / 


( 22 ) 


For distillation, Alice and Bob follow the scheme de¬ 
scribed in Sec. |3.1[ i.e. they start with two qubit pairs 
p = Pf <8> Pf and apply the XY interaction with pulse 
areas a and /3 to their respective qubit pairs. After a 
subsequent measurement of the target qubits, Alice and 
Bob keep the control pair if they obtain equal measure¬ 
ment results. The fidelity F outiXY (F, a, fi) of the source 

3 The square of Vswap 1 is also the SWAP operation and 
it can be understood as another root of SWAP. 


state can be increased depending on the pulse areas a 
and /3, and a formula for F outXY (F, a, p) can be found 
in Ref. HU As discussed before, In the case a = —/?, 
i.e. when both parties apply mutually inverse operations, 
the result coincides with Eq. and the gain in fi¬ 

delity is thus maximal for a = tt/2. In the optimal case 
here, the different qubit interactions correspond to gates 
whose double application result in the isWAP gate. 


3./ Dipole-dipole interaction 


The dipole-dipole coupling between two magnetic mo¬ 
ments pii = ySj and Pj = qS j, separated by a distance 
r, is described by the Hamiltonian [22] 

2 

Hdd = ^3 (S., • Sj — 3 (Si • e r ) (S j ■ e r )). (23) 

Here e r is a unit vector pointing along the connecting 
line between the two identical magnetic moments with 
gyromagnetic ratio 7 and po is the vacuum permeability. 
For example, the electron spins of two nitrogen-vacancy 
centers in diamond that are close enough to each other 
can be coupled via the interaction of the associated mag¬ 
netic moments and entangled in this way [27]. Without 
loss of generality, we can assume the connecting line to 
define the 2 axis and thus, obtain 


H dd = 


P-oT 

107rr 3 


(aWaW+aWaW-toWvW), (24) 


where we assume spin- 1/2 systems that are magneti¬ 
cally coupled. The Hamiltonian H dd is thus obtained 
for anisotropy parameter of 5 = —2. The strength of the 
interaction could in principle be varied by changing the 
distance r between the qubits, which might not be a triv¬ 
ial task. However, as proof of principle of our developed 
concept and to demonstrate that it works for a variety 
of Hamiltonians, we apply the asymmetric distillation 
scheme developed above as well to qubits coupled via 
H dd . We define the pulse area as a = f 0 At' poj 2 /(16nr(t’) 3 ) 
and assume a time-dependent distance r(t). The fidelity 
F out>dd (F,a,/3) after one distillation round with initial 
fidelity F is calculated to be 


-^out,dd(-F, cr, /3) 


v dd (F, a,P) 
S(F,a,p) ’ 


(25) 


and the numerator v dd (F, a, /3) is given by the expression 


^ dd (F, a, P) = (2F(4F + 1) - 1)[2 cos(2(a + p)) 

+ cos(4(a + P)) + 2 cos(6(a + /?))] 

- 2 {F - 1)(4F - 1) cos(4(a - /?)) 

+ 4F(4F+l) + 7. (26) 


Upon detailed inspection, one finds 

-Fout.dd (-F, = F out (F, a, —a), (27) 
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and therefore the discussion of Sec. |3.2| also applies for 
entanglement distillation in case of qubits coupled via 
magnetic dipole-dipole interaction, with the optimal dis¬ 
tillation achieved for a pulse area of a = 7t/8. 


4 Symmetric entanglement distillation with 3 or 
more pairs of spins 

4-1 Extension to three qubit pairs 

In this section we will extend the above setting to a sce¬ 
nario where Alice and Bob have access to three bipar¬ 
tite qubit pairs in a global state p®' 3 and local control 
on isotropic exchange interactions between next nearest 
neighbors, see Fig. [Tj We number the qubits from 1 to 6, 
where Alice has access to odd numbers and can control 
exchange interactions between the qubit pairs (1,3) and 
(3,5). Analogously, Bob has access to even numbered 
qubits and controls interactions between the pairs (2,4) 
and (4,6). 

As in Sec. |3.1[ we will consider protocols where both 
parties first apply controlled sequences of local exchange 
interactions, resulting in overall unitary operations Ua 
and Ub 0 Then they apply a filter based on one round of 
classical communication. This filter is implemented by 
measuring each of the qubits 3 — 6 in the computational 
basis and keeping the state of the qubit pair (1, 2) when¬ 
ever the measurements on the qubit pairs (3,4) and (5,6) 
coincide, see Fig. [3] This is described by the projection 

P,l = (|0 3 0 4 ) (0 3 0 4 | T |1 3 1 4 ) (I 3 I 4 I) 

® (10 5 0 6 > <0 5 0 6 1 + 11 5 1 6 ) (l 5 l e I), (28) 

and the output state p out with fidelity F out relative to the 
maximally entangled target state | T > + ) ( <T + | is obtained 
with success probability P succ , given by 

P succ := tr(I 12 0 Py UaUbpTuIu^) (29) 

Pout := TP— tr 3 , 4)5 , 6 (Ii 2 ® P|| U A U B pfU%U\) (30) 

3UCC 

Pout := TP— tr(|«P+><*+|® P|| U a UbP$ 3 U* b U' a ) . 

* 3UCC ' ' 

(31) 


1 

2 

3 

4 

5 

6 



1 

2 


Figure 3 Circuit of protocols implemented by Alice and 
Bob via iterated applications of controlled exchange interac¬ 
tions on the qubits (1,3,5) and (2,4,6) followed by a mea¬ 
surement on the qubits (3,4,5,6). 


Alice can implement the basic operations, see (18), 


2 i So dt ’ anc J g */o dt' #35(0 


(32) 


by switching on and off an isotropic exchange interaction 
for a specific time t. Up to an irrelevant global phase, 
these operations can be expressed El , as 


U 13 (0 := e *« Fl3 

U 35 (x) :=e^, (33) 


where F,y denotes the flip operation, i.e. it permutes the 
ith and jth tensor factor in (C 2 * )® 6 . If Alice iterates the 
operations (33) with time steps {£,} and {Xi} all uni- 
taries she can implement are of a form 


U A = n tfi3&)tf35(Xi) = n^ Fl3 ^ if35 - (34) 

i i 


4-2 The reachable set of unitaries 

At first we will have to characterize the reachable set 
of unitaries, Ua and Ub, which could be implemented 
by sequences of Alice’s and Bob’s basic operations. This 
characterization can be done separately for the two par¬ 
ties and we will only consider Alice’s side explicitly. 

4 For clarity of notation: the unitaries Ua and Ub are as¬ 

sumed to be represented as matrices on (C 2 )® 6 with Ua act¬ 

ing as identity on Bob’s qubits and Ub as identity on Alice’s. 


The idea for simplifying long products of such op¬ 
erators, with judiciously chosen parameters £ and x is 
to utilize the commutation relations of the flip opera¬ 
tors. Indeed, if the exponentials are expanded, each fac¬ 
tor will be a product of the operators F 43 and F 35 , and 
these can all be evaluated to some permutation oper¬ 
ator of the three sites (1,3,5), i.e., one of the opera¬ 
tors {I,Fi 3 ,F 35 ,Fi 5 ,Zi 3 5 ,Zi 53 }, where Z ijk denotes the 
(anti-)cyclic permutation. These operators span a finite 
dimensional algebra A, for which a convenient basis [2B( 
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is 


A + = - (I + F13 + F35 + F15 + Z135 + Z153) 

A - = - (I — F13 — F35 — F15 + Z135 + Z153) 

Aq = - (21 — Z135 — Z153) 

A \ = - (2F35 — F13 — F15) 

An. = -^= (F13 — F15) 

A3 = (^135 — ^153) • (35) 


Here A + , A_ and Aq are three orthogonal projectors 
summing up to I. They correspond to different irreducible 
representations of the permutation group acting on three 
qubits. A + and A_ are the trivial and the alternating 
representation, which act as projectors on the symmet¬ 
ric and antisymmetric subspace. Aq corresponds to a two 
dimensional representation on which the matrices Ai, A 2 
and A 3 act as Pauli matrices, i.e. [A/,A m ] = 2 iei mn A n 
and A? = A% = A\ = AnTA_i_. 


Now any product of a form as in (34) can be com¬ 


puted in the basis (35) yielding a unitary that is in the 


algebra A. As there is no fully antisymmetric state of 
three qubits we do not further have to take into account 
A_. Hence Aq acts like I — A + such that, up to an irrel¬ 


evant global phase, (34) can always be written as 


U A = e 


_ i(uA + +a-A) 


(36) 


with parameters a £ ( 0 , 27t) and the vectors a £ R 3 such 
as A = (A 1 , A 2 , A 3 ). 

Likewise Bob’s unitaries are described by parameters 
P £ ( 0 , 27t) and bet 3 as 

U B = e i(/3B++bB \ (37) 


with B + and B are defined on the qubits (2,4, 6 ) in the 
same manner as for Alice. 

In our case also the converse holds: Every unitary of 
the form (36) can be obtained as a product as in (34). 
The basic criterion for this is that the operators I, F 13 
and F 35 and their iterated commutators span the whole 
algebra A [291 theorem 2.3]. Finding an explicit and ef¬ 
ficient decomposition is in general a complicated task 
which is the subject of control theory. A good introduc¬ 
tion to this interesting topic can be found in [23130]. 


f. 3 Pretty good protocols 


We can get a qualitative overview over the attainable 
characteristics of possible protocols by random sampling 
of Ua and Us, i.e., by choosing (a,a,/3,b) in (36) at ran¬ 
dom. Fig. [4] shows such a sample of attainable values of 
F out and P succ for different fixed input fidelities F. Good 


protocols in this set are those with a favorable trade-off 
between F out and P s ucc . This can be made precise by the 
notion of Pareto efficiency We say that one proto¬ 
col dominates another whenever it attains higher fidelity 
and a higher success probability, and at least one of these 
parameters is even strictly higher. A protocol which can 
not be dominated by any other is said to be Pareto ef¬ 
ficient, and the corresponding set of pairs (F out , P succ ) 
attained by Pareto efficient protocols is called the Pareto 
front. By definition, the front is a tradeoff curve, along 
which higher fidelity means lower success probability and 
conversely. 



0.872 

0.808 


0.695 

o 

o 

13 

q“ 0.585 
0.538 


0.429 


0.60.625 


0.75 0.797 

Fout 


0.9 0.932 


Figure 4 Random sample of F ou t and P succ for achievable 
protocols with input fidelities F £ (0.6, 0.75, 0.9). The respec¬ 
tive Pareto fronts are marked in red. The Pareto dominated 
region is bounded by dashed red lines. The sample size is 
2 • 10 5 for each case. 


By numerical optimization we identify a family of 
Pareto efficient protocols, as those with parameters 


a = ft = 7 T a = b = (0, 0,r), (38) 

with r £ (0, 7 t/3 ). Remarkably, these protocols are, as in 
the case of two qubit pairs, independent of the input fi¬ 
delity F. The output fidelity and the success probability 
describing the Pareto front is shown in Fig. [5] and can 
be computed as 

16 (4F 2 - 5F + l) cos(3r) + 226 F 2 +F + 16 
out “ 32 (4 F 2 - 5F + 1) cos(3r) + 128F 2 +2F + 113 

(39) 

_ (2F+1) 

succ 729 

x (32 (4F 2 - 5F + 1) cos(3r) + 128E 2 + 2F + 113) . 

(40) 


For r = 7 t /3 the highest success probability and the low¬ 
est output fidelity is attained. In this case the efficiency 
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Figure 5 Output fidelity and success probability of the 
Pareto efficient protocols from ( |40| ) (solid lines) . Output fi¬ 
delity and success probability of the optimal two quit proto¬ 
col (191 iterated on three qubit pairs (dotdashed line). 


equals the case in which Alice and Bob apply no inter¬ 
action at all. In contrast, for r = 0 the highest output 
fidelity and the lowest success probability is attained and 
the maximal achievable fidelity can be computed as 


290F 2 - 79 F + 32 
256F 2 - 158F + 145' 


(41) 


gain with a higher success probability. Moreover also a 
higher fidelity gain is possible when a lower success prob¬ 
ability is accepted. Nevertheless one always has to keep 
in mind that, in the above setting, perfectly controlled 
sequences of interactions are assumed. Hence this might 
harder to realize in an experiment, than an iterated two 
qubit protocol, which can be implemented by only two 
steps of controlled interactions. 


5 Conclusions and Outlook 

We presented entanglement distillation protocols based 
on the exchange interaction using either two or three in¬ 
put pairs. In the case of two input pairs, we analyzed a 
protocol based on (an) isotropic exchange and found that 
entanglement distillation is possible for various interac¬ 
tion types, namely Heisenberg exchange, XY interaction 
and dipole-dipole interaction. If Alice and Bob apply 
mutually inverse operations, it turns out that the output 
fidelity becomes independent of the anisotropy param¬ 
eter £. Further studies could investigate more general 
spin-spin interactions of the form Sf ■ J ■ Sj , with some 
non-diagonal coupling tensor J . An example of such an 
interaction is the so-called Dzyaloshinskii-Moriya inter¬ 
action (321133) . which arises from spin-orbit coupling. 

The above results on three input pairs directly sug¬ 
gest a scheme for finding protocols acting on n-qubit 
pairs by locally controlled next nearest neighbor exchange 
interactions. The operations exp(iaA + ) and exp(i/3B + ) 
can indeed be generalized to arbitrary numbers of qubit 
pairs by choosing A + and B + as projectors on the sym¬ 
metric n-particle subspaces. However, one has to con¬ 
sider that with an increasing number of qubit pairs the 
probability of a joint coincidence of measurements on 
n — 1 qubit pairs decreases exponentially. Hence a more 
detailed investigation is needed to decide whether distil¬ 
lation via exchange interaction can be used to produce, 
with a positive rate, almost maximally entangled pairs 
from a source of sufficiently highly entangled pairs. 
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